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The Anderson transitions in a random magnetic field in three dimensions are investigated numer- 
ically. The critical behavior near the transition point is analyzed in detail by means of the transfer 
matrix method with high accuracy for systems both with and without an additional random scalar 
potential. We find the critical exponent v for the localization length to be 1.45 ± 0.09 with a strong 
random scalar potential. Without it, the exponent is smaller but increases with the system sizes 
and extrapolates to the above value within the error bars. These results support the conventional 
classification of universality classes due to symmetry. Fractal dimensionality of the wave function 
at the critical point is also estimated by the equation-of-motion method. 
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I. INTRODUCTION 

Since the pioneering work by AndersonS, the metal-insulator transition driven by disorder, which is called the 
Anderson transition(AT), has attracted much attention for many yearsoll The critical behavior of the AT is con- 
ventionally classified, depending an the symmetry of hamiltonians, into three universality classes: the orthogonal, the 
unitary and the symplectic classed!! Systems invariant under spin rotation as well as time reversal form the orthogonal 
class. The unitary class is characterized by the absence of the time reversal symmetry. Systems invariant under time 
reversal but having no spin rotation symmetry belong to the symplectic class. 

In the last decade, there has been considerable progress in the numerical study of the AT in three dimensions(3D) 
by the finite-size scaling analysis for quasi-lD systems!! In the early stage, it was not easy to confirm numerically 
for the 3D orthogonal class that the critical exponent is insensitive to the choice of the probability distribution of 
random potential!! This discrepancy in exponents far different distributions of random potential has been,*emoved 
by improving the accuracy of numerical calculationsLi and by taking into account the corrections to scalingcJ. With 
such a high-accuracy analysis, it has been concluded that the critical exponent for the orthogonal system can be 
distinguished from that for the unitary systemllil. These recent developments confirm the universality of critical 
exponents as well as the validity of the conventional classification of universality classes in AT. It should be noted, 
however, that in most cases, such analyses have been restricted to the AT near the band center in the presence of 
a random scalar potential, where the scaling analysis works, fairly well. In contrast, for the AT away from the band 
center, no systematic scaling behavior has been observedcnlj. nn 

The AT in a magnetic field has been studied extensively, mainly in connection with the quantum Hall effectH3'li-i 
Accordingly, in most cases, the magnetic field was assumed to be uniform in space and the disorder was introduced 
by a random scalar potential. On the other hand, in recent years, there has also been considerable interest in the 
transport properties of a system subject to a spatially random magnetic field. The random magnetic field introduces 
randomness as well as the absence of invariance under time reversal in a system. In fact, it has been shownc3 that in 
3D the AT occurs in the presence of the random magnetic field and without a random scalar potential. 

The AT in a random magnetic field is driven by the coherent scattering due to a fluctuating vector potential. A 
nontrivial feature of this coherent scattering by a fluctuating vector potential has been pointed outll3 in a theory of 
strongly correlated spin systems. Much work has also been done on transport properties in 2D in a random magnetic 
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fieldO, in particular in connection with the theory of the fractional quantum Hall effectlia in a high magnetic field. It 
is thus an important issue to understand how the effect of coherent scattering in a strongly fluctuating random vector 
potential will show up in the AT. 

The magnetic field breaks the time reversal symmetry and thus all syste m s i n the magnetic field should belong 
to the unitary class. In fact, it has been demonstrated numerically in 3Dc3E3 that in the presence of a random 
scalar potential, the critical exponent takes a universal value, irrespective of whether the magnetic field is uniform 
or random. The AT with a random potential and in a uniform magnetic field has-been re-analyzed recently and the 
critical exponent for the localization length has been determined to be 1.43 ± 0.06EJ. 

The AT, in 3D, in the presence of a random vector potential a,ud without a random scalar potential, has also been 
investigated based on the finite-size scaling. The data suggested^ that the mobility edge is very close to the band 
edge. The exponent for the localization length has been estimated to be v « lE3 which is considerably smaller than 
that in the case with an additional random scalar potential and in a uniform magnetic field. This seemed to indicate 
that in 3D the AT driven solely by a random vector potential might exhibit critical behavior different from that 
observed in other unitary systems, for example systems having additional random scalar potential. Apparently, this 
questions the validity of the conventional classification of universality classes in AT. On the other iHpd, it should 
be recalled that the finite-size scaling analysis did not work for the AT near the effective band edgecrlla. It is thus 
important to re-examine the applicability of the scaling ansatz to the AT driven solely by the random magnetic field 
in which the mobility edge lies quite close to the band edge. 

In this paper, we report on a high-precision numerical finite-size scaling analysis for the AT in the random magnetic 
field. In order to clarify the origin of the above mentioned discrepancy between the critical exponent of the AT far 
away from the band center induced solely by randomness in a vector potential and the exponent obtained for other 
unitary systems, we have considered systems both with and without an additional random potential. We also evaluate 
the fractal dimension of the wave functions at the critical point based on the equation-of-motion method. 

The paper is organized as follows. In the next section, the hamiltonian which we adopt is introduced. The finite-size 
scaling study on the critical phenomena is presented in section 3. In section 4, the fractal dimensionality of the wave 
function is discussed by means of the equation-of-motion method. Section 5 is devoted to summary and discussion. 



II. MODEL 

The model is defined by the HarniltoniarJll 

H = V eM^CjC, + ]T £iC}Ci, (1) 

<i,j> i 

where C\{Ci) denotes the creation(annihilation) operator of an electron at the site i of a 3D cubic lattice. Energies 
{s^ denote the random scalar potential distributed independently and uniformly in the range [— W/2, W/2]. The 
Peierls phase factors exp(i^j) describe a random vector potential or magnetic field. We confine ourselves to phases 
{Oij} which are distributed independently and uniformly in [— tt,tv]. The hopping amplitude t is assumed to be the 
energy unit, V — 1. The phases {Oij} arc related to the magnetic flux, for example, as 

&id+x + &i+x.i+x+y + Qi+x+y,i+y + &i + V, i = — 27T(fo/</>o, (2) 

where fa and tf>o = hc/\e\ denote the magnetic flux through the plaquette (i, i + x, i + x + y, i + y) and the unit flux, 
respectively. Here x(y) stands for the unit vector in the x(j/)-direction. Note that in the present system, the condition 
that the magnetic flux through a closed surface is zero is satisfied. 



III. FINITE-SIZE SCALING STUDY 

We consider quasi- ID systems with cross section M x AfBi. The Schrodinger equation Htjj = EtJj in such a 
bar-shaped system can be rewritten using transfer matrices T n (2M 2 x 2M 2 ) 

Tpn+l \ _ rr ( i>n \ t _ ( E — H n —I 



(n = 1,2,...) where ip n and H n denote the set of coefficients of the state i\) and the Hamiltonian of the n— th slice, 
respectively. The identity matrix is denoted by /. The off-diagonal parts of the transfer matrix T n can be expressed 
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by the identity matrix because the phases in the transfer-direction can be removed by a gauge transformationE-3. The 
logarithms of the eigenvalues of the limiting matrix T 



t = hmKn^nn^)] 172 " W 

n — >oo 

i=l i=l 

are called the Lyapunov exponents. The smallest Lwapunov exponent Am along the bar is estimated by a technique 
which uses the product of these transfer matricesErQ. The relative accuracies for the smallest Lyapunov exponents 
achieved here is 0.2% for M < 10 and 0.25% ~ 0.3% for M = 12. The localization length £m along the bar is given 
by the inverse of the smallest Lyapunov exponent, £m = 1/Am- 

The assumption of one-parameter scaling for the renormalizcd localization length Am = £m/M implies 

Am = /(£/M), (5) 

where £ = £(E, W) is the relevant length scale in the limit M — > oci. Near the mobility edge E C (W), £ diverges with 
an exponent v as £ ~ x~ v with x = (E — E c )/E c . If the transition is driven by the disorder W at a constant energy, 
x = (W c — W)/W c . At the mobility edge, Am becomes scale-invariant. The quantity Am is a smooth function of E 
and W, and we can expand it as a function of x as 

oc 

A M =A C + Y,MM 1/U x) n . (6) 

71=1 

By fitting our data to the above function, we can determine the critical exponent v and the mobility edge accurately. 
In practice, we truncated the series (0) at the third order(n = 3) and used the standard x 2 -fitting procedures. The 
error bars are estimated by using the Hessian matrix and the confidence interval is chosen to be 95.4%. 

For the the transition at the band center in the presence of a strong random scalar potential, a clear scaling has 
been observed for presently achievable sizes, 6 < M < 12. In fact, all the data (84 points) for M = 6,8, 10, and 12 
in the range 17.8 < W < 19.8 can be successfully fitted by the fitting function (||) up to the 3rd order, which has six 
fitting parameters including the critical point and the critical exponent. We have estimated the critical disorder and 
the exponent v to be W c — 18.80 ±0.04 and v = 1.45±0.09cil The renormalized localization length A c at the critical 
point is 0-5.58 ± 0.003. The error bars of these estimations are at least a factor of 3 smaller than those of the previous 
estimatesEj. 

In contrast, in the absence of the random scalar potential (W = 0) or in the presence of an_additional weak random 
scalar potential (W — 1), for which the critical point lies near the band edge, we have founcO that the correction to 
scaling is not negligible. Near the band edge, the density of states changes rapidly as a function of energy. We have 
thus performed high-accuracy transfer matrix calculations for narrower energy range \E — E c \ < 0.025 around the 
critical point for W = and W = In both cases(W = and W = 1), we have found that the estimation of the 
critical exponent tends to increase with the system-sizes. In order to extrapolate the critical exponent for W = 0, we 
have made calculations for larger system sizes M = 14 and M — 16. Here we show, in Table I, the summary of the 
results for W — obtained by the fittings with different sizes up to M — 16. The relative accuracy in £ M achieved 
for M = 14 and M — 16 is 1% for each sample and 7 and 5 realizations of random phases are considered, respectively. 
The scaling regime is assumed to be [4.39,4.44] as in ref. |24|. It is clear that the critical point exists around E c « 4.41 
(see figure 1). In Table 1, we can see that the exponent v tends to increase with the system-sizes and is likely to 
saturate around v ~ 1.48. Within the error bars, estimated values of v for M > 12 are consistent with 1.45 ± 0.09 
obtained for the band center as well as 1.43 ± 0.06 estimated in the imiform magnetic fieldtil. No evidence has been 
found for cwl which was suggested by calculations with low accuracytJ. The present results support the universality 
of the critical exponent in the unitary systems. The positions of the critical points and the values of A c estimated 
with different combinations of system-sizes are fluctuating for M > 12 (Table I). The value of A c = 0.558 ± 0.003 at 
the band center seems to lie inside the range of this fluctuation. Conventionally, the value of A c is also expected to 
be universal in unitary systems. Our results obtained here seem to be consistent with this universality of A c . 

The mobility edge trajectory in the presence of the random magnetic field is shown in figure 2. Each critical point 
(mobility edge) is estimated based on numerical data by the transfer matrix method with M = 6 ~ 10. It should 
be noted that there exist extended states for energies larger than the critical energy E c w 4.41 for W = 0. This 
type of reentrant phenomena in the ene rgy-d isorder plane has-been commonly observed for systems with the uniform 
distribution of random scalar potentialOcH. It is interpreted^ that the enhancement of extended states for a weak 
additional random scalar potential is due to the enhancement of density of states at that energy regime. 
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IV. EQUATION-OF-MOTION METHOD 



We now turn our attention to the properties of wave function just at-the AT in random magnetic fields. It is well 
known that at the AT, the jKawj -function shows multifractal structurapl which leads to the scale invariant behavior 
of conductance distributionsuJ'EJEil and the energy level statistics.EiTE3 

The direct way to investigate the wave functions is to diagonalize the Hamiltonian. This, however, is numerically 
very intensive. Instead, we calculate here the time evolution of wave packets to extract the information of fractal 
dimension. We first prepare the initial wave packet |0) close to AT by diagonalizing a small cluster located at the 
center of the system. The time evolution of the state at time t is then obtained by 

\t + At) = U(At)\t) 

where U(At) is the time evolution operator. In order to perform effectively the numerical calculation, we approximate 
U(At) by the products of exponential operators as 

U(At) = e - iHAt / h = U 2 ( P At)U 2 ({l - 2p)At)U 2 {pAt) + 0{At 5 ) (7) 
withp= (2- 2 1 / 3 )- 1 and 

U 2 (At) = e - iH i At / 2h . . . e -iir«-iAt/2R e -iff a At/fi e -iH 9 _iAt/2fi _ _ , e ~iH 1 At/2h ^ 

where Hi, ■ ^ H q are decomposition of the original Hamiltonian H = Hi which are simple enough to diagonalize 
analyticallyB3 

The square displacement of the wave packets is defined by 

r 2 {t) = (t\r 2 \t). 

In metallic phase, r 2 {t) is proportional to Dt where D is the diffusion coefficient. In the insulating phase, it saturates 
to the square of localization length, £ 2 . At AT, the anomalous diffusioroo 

r 2 (t) ~ t 2 / d = t 2 / 3 

is expected. The fractal dimension D 2 is estimated from the autocorrelation function 

C{t) = \ f dt'\(t'\0)\ 2 



t _ 

where C(t) is expected to decay as0 

C(t) - r D ^ d . 

In Fig. 3, we show the results of C(t) for the transition at the center of the band in the presence of a strong random 
scalar potentially 7 = 18.8V). By diagonalizing a small cluster of 7 x 7 x 7 located at the the center of the system, 
we follow the time evolution of wave packets in 101 x 101 x 101 systems. Geometric average of C(t) over 10 random 
field and potential configurations are performed. By fitting the data for t > AOTi/V, the fractal dimensionality D 2 is 
estimated to be 

D 2 = 1.52 ±0.18 

considerably smaller than the space dimension d = 3. The above value is consistent with the estimate of 3D system 
at AT in a strong uniform magnetic field .cf] 



V. DISCUSSIONS 



In summary, we have investigated in detail the AT in a random magnetic field based on the transfer matrix method 
with considerably high accuracy. In particular, whether or not the AT driven solely by the random vector potential 
(W = 0) exhibits different critical behavior from other unitary systems has been discussed. In order to clarify the 
above point, we have performed the scaling analysis for the three critical points, namely E — 0, W = 0, and W = 1 
(figure 2). For the transition at the band center (E = 0) in the presence of a strong additional random potential, a 
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clear scaling behavior has been observed and the exponent v has been estimated to be 1.45 ±0.09. This coincides with 
the value obtained for a unitary system in a uniform magnetic fieldlj. It has been found, on the other hand, that the 
correction to scaling is not negligible in the presently achievable sizes for the transitions near the band edge(IU = 1 
and W = 0). The exponents estimated for W — by larger system sizes are consistent with those obtained for other 
unitary systems within the error bars. From the size dependence of v, in contrast to the suggestion in ref. 
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no 

evidence has been found for v « 1. These results indicate the universality of v in the unitary class and hence support 
the conventional classification of the AT by universality classes due to symmetry. 

The mobility edge trajectory has been also obtained in the presence of the random magnetic field. It's qualitative 
shape turns out to be similar to those obtained for other systems with the uniform distribution of random scalar 
potential. 

We have also studied the diffusion of electrons at the AT in the presence of a random magnetic field. By solving the 
time-dependent Schrodinger equation numerically, we examine the time evolution of wave packets at the AT. From 
the asymptotic behavior of the autocorrelation function, we have extracted the fractal dimensionality of the critical 
wave function at the band center. 
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(Mi,M 2 ) Ac v Ec N xL 

(6,8) 0.514 ±0.005 1.05 ± 0.07 4.414 ± 0.001 42 - 84.3 

(8,10) 0.516 ±0.007 1.26 ±0.09 4.414 ± 0.001 42 -30.0 

(10,12) 0.51 ±0.01 1.32 ±0.12 4.414 ± 0.001 42 ~ 15.9 

(12,14) 0.56 ± 0.02 1.49 ±0.16 4.409 ± 0.002 42 -18.7 

(14,16) 0.48 ± 0.02 1.475 ±0.19 4.417 ± 0.002 42 - 32.4 

(12,16) 0.53 ±0.01 1.46 ±0.09 4.413 ± 0.001 42 - 18.6 



TABLE I. Results for A c , v, and E c by the fits using the data for two system-sizes Mi and M2 for W = 0. Here N stands 
for the number of data used for the fit, and Xmin denotes the minimun value of \ 2 obtained by the fit. 
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A 6 




A 

A 8 




Aio 




A12 


A 

A14 




Aie 




4.3900 


0.650 ± 


.001 


0.696 ± 


.001 


n fro/ 1 I r» 

0.736 ± 


.001 


0.777 ± 0.002 


0.808 ± 


.002 


0.843 ± 


.003 


4.3925 


0.634 ± 


.001 


0.675 ± 


.001 


0.709 ± 


.001 


0.745 ± 0.002 


0.768 ± 


.003 


0.799 ± 


.002 


4.3950 


0.620 ± 


.001 


0.654 ± 


.001 


0.686 ± 


.001 


0.713 ± 0.002 


r\ 'ran i t~\ 

0.729 ± 


.002 


0.765 ± 


.003 


4.3975 


0.605 ± 


.001 


0.634 ± 


.001 


0.658 ± 


.001 


0.683 ± 0.002 


0.696 ± 


.002 


0.721 ± 


.004 


4.4000 


0.588 ± 


.001 


0.614 ± 


.001 


0.636 ± 


.001 


0.654 ± 0.002 


0.664 ± 


.002 


0.688 ± 


.002 


4.4025 


0.574 ± 


.001 


0.595 ± 


.001 


0.609 ± 


.001 


0.627 ± 0.002 


0.635 ± 


.002 


0.652 ± 


.005 


4.4050 


0.561 ± 


.001 


0.575 ± 


.001 


0.588 ± 


.001 


0.600 ± 0.002 


0.603 ± 


.001 


0.613 ± 


.004 


4.4075 


0.546 ± 


.001 


0.557 ± 


.001 


0.565 ± 


.001 


n rwr i r\ r\r\-t 

0.575 ± 0.001 


0.573 ± 


.002 


0.587 ± 


.001 


4.4100 


0.530 ± 


.001 


0.540 ± 


.001 


0.546 ± 


.001 


C\ r— A j \ 1 / \ Al Al 1 

0.549 ± 0.001 


0.547 ± 


.002 


0.554 ± 


.003 


4.4125 


0.517 ± 


.001 


0.522 ± 


.001 


0.524 ± 


.001 


0.526 ± 0.001 


0.522 ± 


.002 


0.529 ± 


.002 


4.4150 


0.506 ± 


.001 


0.505 ± 


.001 


0.505 ± 


.001 


0.502 ± 0.001 


0.498 ± 


.002 


0.497 ± 


.002 


4.4175 


0.494 ± 


.001 


0.492 ± 


.001 


f\ A Cl A 1 l\ 

0.484 ± 


.001 


r\ a o "1 1 / \ /~1A1 1 

0.481 ± 0.001 


0.474 ± 


.002 


0.473 ± 


.001 


A A C\f\f\ 

4.4200 


0.483 ± 


.001 


f\ A A i r\ 

0.474 ± 


.001 


0.468 ± 


.001 


0.461 ± 0.001 


0.453 ± 


.001 


0.452 ± 


.002 


4.4225 


0.472 ± 


.001 


0.459 ± 


.001 


0.450 ± 


.001 


o ^ j o i r\ r\r\-\ 

0.443 ± 0.001 


0.432 ± 


.001 


f\ A C\ A 1 ("1 

0.424 ± 


.002 


4.4250 


0.457 ± 


.001 


r\ a a t~ i t~\ 

0.445 ± 


.001 


0.433 ± 


.001 


C\ A c\<~\ i c\ r\r\-t 

0.423 ± 0.001 


0.411 ± 


.001 


0.404 ± 


.002 


4.4275 


0.448 ± 


.001 


0.431 ± 


.001 


0.417 ± 


.001 


0.405 ± 0.001 


0.393 ± 


.002 


0.383 ± 


.001 


4.4300 


0.438 ± 


.001 


0.418 ±0 


.001 


0.401 ± 


.001 


0.387 ± 0.001 


0.373 ± 


.001 


0.364 ± 


.001 


4.4325 


0.428 ± 


.001 


0.405 ± 


.001 


0.387 ± 


.001 


0.371 ± 0.001 


0.355 ± 


.001 


0.345 ± 


.001 


4.4350 


0.420 ± 


.001 


0.391 ±0 


.001 


0.374 ± 


.001 


0.356 ±0.001 


0.339 ± 


.001 


0.328 ± 


.001 


4.4375 


0.409 ± 


.001 


0.382 ± 


.001 


0.359 ± 


.001 


0.342 ± 0.001 


0.323 ± 


.001 


0.313 ±0 


.001 


4.4400 


0.402 ± 


.001 


0.369 ± 


.001 


0.347 ± 


.001 


0.328 ±0.001 


0.308 ± 


.002 


0.298 ± 


.001 



TABLE II. List of numerical data of Am presented in figure 1 for W — 0. 
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0.9 
0.8 
0.7 

0.6 
0.5 

0.4 
0.3 



4.39 4.4 4.41 4.42 4.43 4.44 

E 

FIG. 1. Renormalized localization lengths for W = are presented in log-scale as a function of energy. The dots, the crosses, 
the triangles, the squares, the diamonds, the circles correspond to M = 6, 8, 10, 12, 14, and 16, respectively. 




FIG. 2. Mobility edge trajectory for the 3D system in the random magnetic field. The re-entrant phaenomenon is clearly 
seen. 
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